Gapless two-flavor color superconductor 



Igor Shovkovy^'* and Mei Huang^'^ 
^ Institut fiir Theoretische Physik, J.W. Goethe-Universitdt, D-60054 Frankurt/Main, Germany 
^Physics Department, Tsinghua University, Beijing 100084, China 
(February 1, 2008) 



cn 
o 
o 

(N 

Oh- 
< 

(N 
m ■ 

> : 

^ ■ 

O 

m 

o : 

i!^ ■ 
Oh: 

Oh; 
^ ■ 



X 



A new, gapless two-flavor color superconducting phase that 
appears under conditions of local charge neutrality and (3- 
equilibrium is revealed. In this phase, the symmetry of the 
ground state is the same as in the conventional two-flavor 
color superconductor. In the low-energy spectrum of this 
phase, however, there are only two gapped fermionic quasi- 
particles, and the other four quasiparticles are gapless. The 
origin and the basic properties of the gapless two-flavor color 
superconductor are discussed. This phase is a natural candi- 
date for quark matter in cores of compact stars. 

PACS numbers: 12.38.-t, 26.60.+C 



Sufficiently cold and dense quark matter is a color su- 
perconductor [1,2]. The ground state is characterized by 
a condensate of Cooper pairs that are made of quarks 
with equal and opposite momenta. In this phase of mat- 
ter, the SU(3)c color gauge group is broken (partially 
or completely) through the Anderson-Higgs mechanism. 
At asymptotic densities, this phenomenon was studied in 
detail from first principles in Refs. [3-5]. 

In general, quark matter at high baryon density has 
a rich phase structure, consisting of different normal 
and color-superconducting phases. Changing the baryon 
chemical potential, it should be possible to obtain phases 
in which either two or three of the lightest quark flavors 
participate in Cooper pairing. 

At not very large densities, one gets quark matter made 
of only up and down quarks. If the Fermi momenta of 
different types of quarks are approximately the same, 
the up and down quarks will form Cooper pairs in the 
color-antitriplet, flavor-singlet, spin-zero channel [2-4]. 
The corresponding ground state of matter is the so-called 
two- flavor color superconductor (2SC). In this phase, the 
color gauge group is broken by the Anderson-Higgs mech- 
anism down to SU(2)c subgroup. With the conventional 
choice of the condensate pointing in the "blue" direction, 
one finds that the condensate consists of red up (ur) and 
green down (dg), as well as green up (ug) and red down 
(dr) diquark Cooper pairs. The other two quarks (ub and 
db) do not participate in pairing. This is the conventional 
picture of the 2SC phase [2-4]. 

At sufficiently large densities, i.e., when the value of 
the chemical potential exceeds the constituent (medium 
modified) mass of the strange quark, quark matter should 
consist of all three quark flavors. This opens the possibil- 
ity for Cooper pairing between the three lightest quarks 
(up, down and strange). Depending on the value of 



the strange quark mass, as well as other parameters in 
the theory, one might get the color-flavor locked phase 
(CFL) [6], or even the more exotic crystalline phase [7] 
of quark matter. Besides that, it is also possible that 
up and down quarks form 2SC matter, while the strange 
quarks do not participate in pairing. We call this latter 
phase 2SC+S. Still another possibility is that the strange 
quarks by themselves form a color superconducting con- 
densate. This will be then a spin-1 condensate [4,8]. 

It is natural to expect that some color superconducting 
phases may exist in the interior of compact stars. The 
estimated central densities of such stars might be as large 
as lOpo (where po « 0.15 fm~'^ is the saturation density), 
while their temperatures should be in the range of tens 
of keV. These values are very encouraging. 

Matter in the bulk of compact stars should be neutral 
with respect to electric as well as color charges. Also, 
such matter should remain in /3-equilibrium. Satisfying 
these requirements may impose nontrivial relations be- 
tween the chemical potentials of different quarks. In turn, 
such relations could substantially influence the pairing 
dynamics between quarks, for instance, by suppressing 
some color superconducting phases and by favoring oth- 
ers. 

It was argued in Ref. [9], that the 2SC+S phase be- 
comes less favorable than the CFL phase if the charge 
neutrality condition is enforced (note that the strange 
quark mass was chosen too small to allow the appearance 
of a pure 2SC phase in Ref. [9]). Therefore, one might 
speculate that only CFL quark matter could exist inside 
compact stars. In such a case, the outside hadronic layer 
of the star would make a direct contact with the CFL 
quark core through a separating sharp interface [10,11]. 
The baryonic density and the energy density would have 
a jump (smoothed only over microscopic distances) at 
the interface. Of course, a mixed phase of hadronic and 
CFL matter is another possibility, but it does not seem 
to be favorable [11] (for other studies of the stars with 
CFL quark matter in their interior see Ref. [12]). 

We should mention that the general conclusion of Ref. 
[9] was essentially confirmed in Ref. [13] where it was 
claimed that the 2SC-I-S could exist only in a narrow 
window (about 10 to 15 MeV wide) of baryon chemical 
potential around the midpoint /i = /^s/3 « 450 MeV. For 
lower values of fi, no neutral 2SC phase was found. In 
this Letter, however, we argue that a neutral 2SC phase 
does exist at lower values of fi. It is a gapless rather than 
an ordinary 2SC phase. In particular, the conventional 
density relations between the pairing quarks, used in Ref. 
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[13], are not valid in this new type of 2SC quark matter. 
This phase is rather unusual because it has the same 
symmetry of the ground state as the conventional 2SC 
phase, but the spectrum of its fermionic quasiparticles is 
very different. 

It appears that the neutral gapless 2SC phase of quark 
matter has already been used in numerical studies of 
Refs. [14,15] (and, possibly, in Ref. [16]). Its special 
properties, however, have not been appreciated and have 
never been discussed before. In this Letter, we are going 
to explain the origin of this new phase and shed some 
light on its properties. 

Let us start our analysis from discussing the quark 
model that we use. Without loosing generality, we as- 
sume that the strange quark is sufficiently heavy and 
does not appear at intermediate baryon densities under 
consideration (e.g., this might correspond to the baryon 
chemical potential = /Ub/S in the range between about 
350 and 450 MeV). Then, in our study, we could use the 
simplest SU(2) Nambu Jona-Lasinio (NJL) model [14]. 
The explicit form of the Lagrangian density reads: 

C = q{i^^'^^ - mo)q + Gs [{qqf + {qijsrqf] 

+Gn [(*g^£e''75a)(ig«''759^)] , (1) 

where q'~' = Cq^ is the chargc-conjiigate spinor and 
C = «7^7° is the charge conjugation matrix. The quark 
field q = qia is a four-component Dirac spinor that car- 
ries flavor (i = 1,2) and color (a = 1,2,3) indices, 
r = (t^,t^,t'^) are Pauli matrices in the flavor space, 
while (e)*'^ = e*'^ and {e'')°'^ = e"-^^ are antisymmetric 
tensors in flavor and color, respectively. We also intro- 
duce a momentum cut-oflP A, and two independent cou- 
pling constants in the scalar quark-antiquark and scalar 
diquark channels, Gs and Gd- 

The values of the parameters in the NJL model are 
chosen as follows; Gs = 5.0163 GeV-^ and A = 0.6533 
GeV [14]. In this Letter, we consider only the chiral 
limit with mo = 0. As for the strength of the diquark 
coupling G/j, its value is taken to be proportional to the 
quark-antiquark coupling constant, i.e., Gd — vGs with 
a typical number for rj being around 0.75 [13,14]. 

In /3-equilibrium, the diagonal matrix of quark chemi- 
cal potentials is given in terms of baryonic, electric and 
color chemical potentials. 



Hub = M - g/^e 



(5) 

(6) 



The effective potential for quark matter at zero tem- 
perature and in /3-equilibrium with electrons takes the 
form 
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where f2o is a constant added to make the pressure of 
the vacuum zero. This is the zero temperature limit of 
the potential derived in Ref. [14]. The sum in Eq. (7) 
runs over all (6 quark and 6 antiquark) quasiparticles. 
The dispersion relations and the degeneracy factors of 
the quasiparticles read 



E^, = E{p) ± 
iJ± = E(p) ± ^lab, 



V[S(p)±Ai]' + A2±^/i, 



[xl] 
[xl] 
[x2] 
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where the shorthand notations are E(j)) = y'p^ -|- w?, 

p, = {Hur + Hdg)/'^ = A* - A«e/6 + /ig/S, and 5lX = {jldg - 

/^„r)/2 = A*e/2. 

In this study we concentrate exclusively on the 2SC 
phase of quark matter. The numerical study shows that 
the constituent quark mass m is zero in this phase [14]. 
Then by making use of the dispersion relations, we derive 
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127r2 AGd 27r2 127r2 127r2 
-2^ ^ + Ji? + A2 + ^{p - n? + A2) 

-2e{Sf,-A)J^_ ^(<5m-V(p-m)2 + A2), (11) 

where = /l ± -y/ ((5/i)^ — A^. Note that the physical 
thermodynamic potential that determines the pressure, 
f^phys = —P, is obtained from O in Eq. (11) after substi- 
tuting /ig , /ie and A that solve the color and electrical 
charge neutrality conditions, as well as the gap equation, 
i.e., 



l^ij,a/3 = ifiSij - IJ-eQij)Sa/3 + -^I^S^ij{T8)a/3, (2) 

where Q and Tg are generators of U(l)em of electromag- 
netism and the U(l)8 subgroup of the color gauge group. 
The explicit expressions for the quark chemical potentials 
read 



_ _ 2 1 

Mur — fJ-ug — M ~ gA'e + gA^S) 

1 1 

Mdr = Mdg = M + g/^e + 3M8, 



(3) 
(4) 



dn dn dn , , 

«8 = — =0, nQ = — =0, and ^=0. (12) 



Let us assume that the solution to the color neutrality 
condition, /X8(/Lie,A), is known for a given value of fj.. 
In practice, we construct this function using a numerical 
set of solutions. By substituting this solution into the 
potential in Eq. (11), one is left with a function Q that 
depends only on two parameters, He and A. In this two- 
dimensional parameter space, the solution to each of the 
extra two conditions in Eq. (12) is a one-dimensional 
line. Graphically, this is shown in Fig. 1. Since both 
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conditions should be satisfied simultaneously, one needs 
to find the intersection points of the corresponding lines. 

First we discuss the neutrality condition, ng = 0. The 
corresponding line of solutions in the two-dimensional pa- 
rameter space is a monotonic, single valued function. In 
Fig. 1, it is represented by the thick dash-dotted line. 
This line of solutions consists of two branches in two qual- 
itatively different regions of the parameter space that are 
separated by the line A = In the leading order ap- 

proximation, this neutrality line is independent of the 
coupling constant Go- This should be obvious from the 
derivation of the first two conditions in Eq. (12). 



A ( MeV) 




FIG. 1. The graphical representation of the solutions to the 
gap equation for three different values of the diquark coupling 
constant (thick solid and dashed lines), and to the electrical 
neutrality condition (thick dash-dotted line) . The intersection 
points represent the solutions to both equations. The thin 
solid line divides two qualitatively different regions, A < /.te/2 
and A > /ie/2. The results are plotted for ^ — 400 MeV and 
three values of diquark coupling constant Gd = V^s with 
77 = 0.5, 77 = 0.75, and tj = 1.0. 

The line of solutions to the gap equation is most inter- 
esting. It is made of two different branches in two regions, 
A < fj,e/2 and A > /ie/2, see Fig. 1. The upper is the 
main branch that exists down to /ig = (i.e., no mis- 
match in the Fermi surfaces of the up and down quarks). 
The lower branch appears only in a finite window of elec- 
tric chemical potentials, and it merges with the upper one 
at a point on the line A = fie/2. In Fig. 1, we show three 
solutions to the gap equation a,t fi = 400 MeV in three 
regimes with different coupling constants, corresponding 
to 77 = 0.5, r] = 0.75, and r] = 1.0. [The results look the 
same at all values of The three curves have qualita- 
tively the same shape, but differ by overall scale factors. 
The difference in the overall scale has an important con- 
sequence. At weak coupling (r/ ^ 0.7), there is no neutral 
2SC phase because there is no intersection of the solution 
to the gap equation with the neutrality line, see t] — 0.5 
solution in Fig. 1. At intermediate (0.7 ^ 77 < 0.8) and 
strong (?7 > 0.8) coupling, on the other hand, there is 
an intersection with the neutrality line at a point on the 
lower and upper branches, respectively. In these last two 
cases, neutral 2SC phases exist. 



If the electric chemical potential were a free parameter, 
one would find that increasing its value leads to a first 
order phase transition [17]. This is the result of the ap- 
pearance of two competing local minima (at A = and 
A 7^ 0) of the potential r2(A, — const) for a range of 
intermediate values of the parameter /Xg. Note, that the 
location of the local maximum of the potential that lies 
between the two local minima is determined by the lower 
"unstable" branch of the solutions to the gap equation 
in Fig. 1. The first order phase transition happens at 
about fif = \/2Ao where Aq is the value of the gap at 
fie — (this estimate is derived in the approximation of 
weak coupling [18]). From our discussion below, it will be 
clear that this first order phase transition is unphysical 
under the requirement of local neutrality of quark mat- 
ter. This is because the transition typically happens be- 
tween two types of matter with different, nonzero charge 
densities (e.g., positively charged color superconducting 
matter and negatively charged normal quark matter). It 
is worth mentioning, however, that this first order phase 
transition can get physical meaning in a globally neutral 
mixed phase of quark matter [19], and in some condensed 
matter systems where there is no analogue of the charge 
neutrality condition [20]. 

The three-dimensional view of the potential as a 
function of two parameters, A and /ig, is represented by 
a surface in Fig. 2. The results are plotted for /i — 400 
MeV and for the diquark coupling with 77 = 0.75. In 
this study, we are interested exclusively in locally neutral 
quark matter (mixed phases are discussed in Ref. [19]), 
and therefore we consider the potential only along the 
neutrality line (represented by a black solid line in Fig. 
2). 




FIG. 2. The effective potential as a function of the diquark 
gap A and electric chemical potential /j,^. The black solid 
line gives the potential along the electric neutrality line. The 
white solid line shows the line of solutions to the gap equation. 
The results are plotted for /i — 400 MeV and Gd ~ ^Gs with 

= 0.75. 

The study shows that, for any coupling strength, such 
a potential of the neutral quark matter, Q{A, fie\nQ=o), 
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as a function of A has one local minimum which is also 
the global one. In the regime of intermediate strength 
of diquark coupling (e.g., i] = 0.75), this minimum cor- 
responds to a point on the lower branch of the solutions 
to the gap equation in Fig. 1 (i.e., A « 68 MeV and 
fie ~ 148.4 MeV). This may sound very surprizing be- 
cause the points on the lower branch correspond to local 
maxima of the potential V{A, He = const). To make this 
clear, we compare the effective potential calculated for 
a fixed value of jJe with the potential defined along the 
neutrality line in Fig. 3. As we see, after enforcing the 
condition of charge neutrality, a local maximum of the 
former becomes a global minimum of the latter. 

It is interesting to study the dependence of the diquark 
coupling constant on the properties of the potential of 
neutral matter. In the weakly coupled regime (77 < 0.7), 
the potential has the global minimum at A = 0. This is in 
agreement with the fact that the neutrality line does not 
intersect the corresponding line of nontrivial solutions 
(e.g., 77 = 0.5 case in Fig. 1). For intermediate values of 
the coupling (0.7 ^ rj ^ 0.8), the extremum at A = is 
a local maximum, while a minimum appears at A ^ 0, 
see solid curve in Fig. 3. In this case, the minimum 
corresponds to an intersection of the neutrality line with 
the lower branch of the solutions to the gap equation 
(e.g., 77 = 0.75 case in Fig. 1). Finally, in the strong 
coupling regime (77 > 0.8), there is a nontrivial minimum 
which corresponds to an intersection of the neutrality 
line with the upper branch of the solutions (e.g., r] = 1.0 
case in Fig. 1). The phase transition controlled by the 
coupling constant 77 is a second order phase transition. 



two gapless and four gapped quasiparticles. The num- 
ber of gapless modes changes when one crosses the line 
A ^ The quasiparticle dispersion relations, orig- 

inating from the red and green quarks, as described by 
Eq. (10), are shown in Fig. 4. Each line corresponds to 
two degencirate quasiparticles. The dispersion relations 
of gapless blue quarks, given in Eqs. (8) and (9), are not 
shown in Fig. 4. 

In ordinary 2SC matter (with a solution in the region 
A > /Xe/2), there are two doublets of gapped modes: 
one with the gap A — and the other with the gap 
A -I- /ie/2. As one approaches the boimdary between the 
regions, A Me/2, the two quasiparticles of the first 
doublet with the smaller gap gradually become gapless. 
They also remain gapless in the phase with A < /Xe/2 
(regime of intermediate coupling). This justifies the 
name "gapless 2SC phase" . 

In the study of the CFL phase in Ref. [21], it was re- 
vealed that the number densities of the pairing quarks 
are equal. By applying the same arguments in the case 
of 2SC quark matter, one might expect to get similar 
relations for the number densities of red up and green 
down (as well as green up and red down) quarks. More- 
over, these relations have been numerically confirmed for 
the 2SC phase of quark matter in the strongly coupled 
regime in Ref. [22]. We find, however, that the num- 
ber densities of pairing quarks are not always equal. In 
particular, they are not equal in the gapless 2SC phase, 
i.e., Uur = nug 7^ Udr = ndg- This is directly related to 
the fact that the corresponding solution is found in the 
region A < /ie/2, where two extra gapless modes appear. 
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FIG. 3. The effective potential as a function of tlie diquark 
gap A calculated at a fixed value of the electric chemical po- 
tential /Xe = 148.4 MeV (dashed line), and the effective poten- 
tial defined along the neutrality lino (solid line). The results 
are plotted for fi = 400 MeV and Go = r?Gs with q = 0.75. 

Here we argue that the regime of intermediate cou- 
plings describes a new, gapless phase of 2SC quark mat- 
ter. This phase of superconducting quark matter pos- 
sesses four gapless and only two gapped fermionic quasi- 
particles in its spectrum, see Eqs. (8)-(10). Recall that 
the ordinary 2SC phase (which, in fact, appears in the 
region A > /Xe/2 in the strong coupling regime) has only 
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FIG. 4. Three types of dispersion relations of quasiparti- 
cles that originate from quarks of two (red and green) colors 
participating in Cooper pairing. The two solid (dashed) lines 
represent the dispersion relations in the regime with A > /Lie/2 
(A < /Je/2). The dash-dotted hne corresponds to the special 
case with A = Each line corresponds to two degenerate 

quasiparticles. 

In conclusion, the analysis of this Letter confirms that 
the charge neutrality condition can strongly affect color 
superconducting quark matter, favoring some phases and 
disfavoring others. Moreover, as we argued, by impos- 
ing the charge neutrality condition, one can also obtain 



4 



new stable phases of matter which could not exist other- 
wise. The gapless 2SC phase of quark matter is a beau- 
tiful example of such a new phase. The symmetry of its 
ground state is the same as that of the conventional 2SC 
phase. However, the energy spectrum of the fermionic 
quasiparticles in gapless 2SC matter is different: it has 
two additional gapless modes. We should note that this 
phase is very different from another example of gapless 
color superconductivity, discussed in Ref. [24], which is 
a metastable state of CFL quark matter. The gapless 
2SC phase is a stable, neutral state of two-flavor quark 
matter. 

In nature, gapless 2SC quark matter could exist in 
compact stars. Indeed, this phase is neutral with re- 
spect to electric and color charges and satisfies the (3- 
equilibrium condition by construction. Of course, satisfy- 
ing these requirements is necessar\^, but not siiflicicnt. In 
order to decide whether the gapless 2SC phase is likely to 
appear inside stars, further detailed studies are needed. 

We could speculate that the thermodynamic properties 
of the gapless 2SC phase should be closer to the proper- 
ties of neutral normal quark matter rather than to those 
of strongly coupled 2SC matter. Qualitatively, this is 
expected from observing the relative shallowness of the 
effective potential in Fig. 3 defined along the neutrality 
line. It indicates that, for a given value of the baryon 
chemical potential /x, the pressure difference of (neutral) 
normal and gapless 2SC quark matter is substantially 
smaller than the pressure difference of the normal quark 
and conventional 2SC phase, 5P = (/iA/Tr)^ [21,23]. 

Here we studied only the most general properties of the 
gapless 2SC phase of quark matter. By taking into ac- 
count that this phase is a realistic candidate for matter 
inside compact star cores, it would be very interesting 
to study its other properties that could potentially af- 
fect some compact star observables. In particular, this 
includes magnetic properties, neutrino emissivities, and 
various transport properties of gapless 2SC matter. 
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